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1  Introduction 


Background 

Increasing  computer  speed  and  the  development  of  new  mathematical 
representations  of  physical  processes  have  outdated  many  water  quality 
models.  Several  models  used  in  the  past  to  simulate  nonpoint  source 
(NPS)  pollution  are  adequate,  but  many  areas  need  to  be  upgraded  to  take 
advantage  of  recent  developments  in  modeling  and  computer  speed. 

Models  are  used  to  help  in  long-term  planning  to  determine  the  best 
methods  to  reduce  NPS  pollution  and  maintain  water  quality.  Reduction  of 
NPS  pollution  is  accomplished  by  using  best  management  practices 
(BMP’s).  Accurate  modeling  is  essential  in  determining  how  various 
BMP’s  will  impact  water  quality,  alone  or  in  conjunction  with  others.  The 
purpose  of  this  report  is  to  compare  and  evaluate  existing  hydrologic  and 
watershed  water  quality  models.  Recommendations  for  a  base  model(s)  to 
be  used  to  predict  NPS  pollution  will  be  made.  Areas  for  model  improve¬ 
ment  will  be  identified  to  increase  the  model’s  reliability  and  usability. 
Viable  modeling  approaches  today  typically  are  comprehensive.  These 
models  contain  components  to  account  for  both  NPS  and  point  source 
pollutants. 

Since  NPS  pollution  is  most  often  associated  with  runoff  events,  the 
hydrologic  portion  of  each  model  will  be  examined  first.  Many  early  mod¬ 
els  used  empirical  relationships  to  determine  runoff  and  peak  flow,  such 
as  the  Soil  Conservation  Service  (SCS)  curve  number,  rational  method, 
SCS  TR-55.  Increases  in  computational  speed  have  allowed  many  models 
to  be  developed  that  represent  each  process  of  the  hydrologic  cycle  indi¬ 
vidually  with  a  model  based  on  the  physics  of  the  process.  The  differing 
manners  in  which  the  models  simulate  infiltration,  evapotranspiration,  per¬ 
colation,  and  runoff  routing  will  be  examined.  The  applicability  of  the 
model  in  both  urban  and  rural  settings  will  be  addressed. 

Initially,  sediment  was  the  only  component  of  NPS  pollution  studied. 
Many  models  use  the  Universal  Soil  Loss  Equation  (USLE)  or  some  modi¬ 
fication  to  predict  sediment  loading.  The  state  of  the  art  now  includes  model¬ 
ing  the  processes  of  detachment,  deposition,  and  transport  rather  than 
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using  the  regression  approach  of  the  USLE.  This  allows  the  user  to 
examine  the  temporal  as  well  as  spatial  changes  in  sediment  loading. 
Increasing  concerns  of  pesticides,  nutrients,  and  other  chemicals  entering 
the  water  system  have  caused  their  inclusion  in  many  watershed  models. 
When  these  other  pollutants  are  considered,  adsorption,  degradation,  and 
solubility  must  also  be  considered.  The  other  pollutants  should  be 
accounted  for  from  the  time  they  enter  the  watershed  until  they  leave  or 
are  degraded,  including  in  channels  and  impoundments.  The  effects  of 
these  pollutants  on  groundwater  should  also  be  taken  into  account. 

Two  items  affect  model  accuracy  more  than  any  others,  scale  and  time 
step.  When  models  are  applied  over  too  large  an  area,  many  parameters 
are  lumped  or  averaged  for  the  area.  This  neglects  the  small-scale  hetero¬ 
geneities  that  are  found  in  the  real  world  and  affect  runoff  quality  and 
quantity.  If  the  time  step  is  too  long,  many  physically  based  models  do  not 
operate  properly,  and  the  benefits  of  using  distributed  rather  than  lumped 
parameters  are  lost.  The  ideal  solution  would  be  to  use  a  distributed 
parameter  model  with  a  short  time  step,  but  this  introduces  two  more  con¬ 
cerns,  data  collection  and  computational  time.  Collection  of  the  necessary 
meteorologic  data  and  watershed  parameters  for  a  long-term  simulation  of 
a  large  watershed  could  be  burdensome  and  quite  voluminous.  To  help 
alleviate  this  burden,  linking  the  model  with  a  geographic  information  sys¬ 
tem  (GIS)  and  using  weather  simulators  is  desirable.  To  reduce  computa¬ 
tional  time,  changes  in  the  time  step  may  be  used  during  times  with  no 
rainfall  to  increase  the  efficiency  of  the  model. 


Water  Quality  Models 

Models  vary  in  many  ways:  time  step,  scale,  whether  the  model  simu¬ 
lates  single  events  or  on  a  continuous  basis,  and  how  different  compo¬ 
nents  are  computed.  For  NPS  modeling,  the  only  feasible  option  is  to 
incorporate  a  continuous  approach.  Loadings  from  a  watershed  area  need 
to  be  represented  over  time,  not  just  for  a  single  event.  The  first  area  that 
will  be  examined  in  this  model  review  is  how  runoff  is  computed.  One 
common  method  to  compute  the  amount  of  runoff  from  daily  rainfall  data 
is  to  use  the  SCS  curve  number  method.  This  empirical  method  is  based 
on  numerous  tests  relating  soil  moisture,  soil  type,  and  land  use  to  the 
depth  of  runoff  from  a  given  depth  of  rainfall.  Other  methods  require  the 
modeling  of  rainfall  abstractions,  such  as  interception,  infiltration,  and 
evaporation,  to  determine  the  amount  of  runoff.  Various  methods  of  differ¬ 
ing  complexity  can  be  used  to  compute  the  abstraction.  Some  models  keep 
track  of  water  once  it  has  percolated  into  the  soil  to  examine  interflow 
and  groundwater  recharge.  Groundwater  percolation  is  a  necessary  compo¬ 
nent  for  most  reasonably  sized  watersheds. 

Once  the  amount  of  runoff  has  been  computed,  routing  of  the  runoff  to 
a  receiving  body  becomes  the  focus.  Some  models  offer  accurate  but  com¬ 
plex  methods  of  solving  the  continuity  and  conservation  of  momentum 
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equations.  These  models  incorporate  either  modified  forms  of  St.  Venant’s 
equation  or  an  approach  known  as  the  kinematic  wave  approximation. 
These  methods  require  a  shorter  time  step  than  is  practical  for  some  mod¬ 
els.  Use  of  synthetic  hydrographs  based  on  storm  duration  and  amount  of 
runoff  may  be  used.  Several  methods  have  been  developed  by  the  SCS  to 
determine  peak  runoff  rate  (TR-20,  TR-55,  and  rational  method).  Many 
routing  methods  are  currently  available  including  those  which  are  fully 
hydrodynamic.  After  the  runoff  has  reached  a  channel  or  reservoir,  it  must 
be  routed  to  the  outlet  of  the  watershed.  Many  routing  methods  are  avail¬ 
able,  including  those  that  are  fully  hydrodynamic. 

Sediment  yield  can  be  computed  in  several  ways.  Most  models  use  a 
variation  of  the  USLE.  The  original  USLE  was  developed  to  predict  long¬ 
term  soil  loss,  but  modifications  allow  its  use  for  individual  storms.  The 
main  variation  between  models  is  whether  rainfall,  runoff,  or  a  combina¬ 
tion  provide  the  energy  for  erosion.  In  urban  models  simple  build-up/ 
wash-off  relationships  are  utilized  to  approximate  erosion.  More  detailed 
models  compute  detachment  based  on  both  rainfall  and  runoff  energy  and 
then  compute  transport  capacity  based  on  runoff  to  determine  sediment 
yield.  Models  may  compute  sediment  deposition  and  erosion  in  channels 
and  ponds. 

If  chemical  transport  is  considered,  several  items  need  to  be  reviewed. 
For  nutrients,  nitrogen  and  phosphorus,  the  biochemical  processes  that 
affect  the  form  of  the  nutrient  need  to  be  taken  into  account.  For  pesti¬ 
cides,  decay  should  be  included  as  well  as  simulation  of  decay  compo¬ 
nents.  Dividing  the  chemical  transport  between  adsorbed  and  soluble 
portion  must  be  considered.  Chemical  effects  on  groundwater  may  be 
included.  Biochemical  processes  need  to  be  simulated  during  all  phases  of 
transport,  not  just  on  the  field.  Crop  growth  and  time  of  year  (impacts 
through  other  time-dependent  processes)  will  affect  chemical  transport; 
therefore,  their  impacts  need  to  be  simulated  as  well. 


Objective 

The  objective  of  this  study  is  to  review  existing  watershed  water  qual¬ 
ity  models  and  select  a  building  block  for  the  development  of  a  Corps  of 
Engineers  modeling  capability  for  predicting  watershed  runoff  water  qual¬ 
ity  as  influenced  by  land  use  change,  development,  watershed  improve¬ 
ment  strategies,  and  BMP’s.  Other  factors  will  be  examined  during  the 
review  of  existing  models.  The  methods  used  to  assign  various  parameters 
to  different  locations  in  the  watershed  will  be  examined.  Input  require¬ 
ments  and  possible  simulation  of  inputs  will  also  be  reviewed. 
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2  Model  Review 


The  following  is  a  brief  overview  of  some  of  the  existing  watershed 
water  quality  models  that  are  readily  available.  The  review  is  in  no  way 
all  inclusive  but  represents  models  that  are  commonly  used  and  well  vali¬ 
dated.  The  water  quality  models  may  be  broken  into  two  groups,  urban 
and  nonurban.  The  urban  models  are  STORM,  SWMM,  and  DR3M- 
QUAL.  Nonurban  models  include  CREAMS/GLEAMS,  EPIC,  SWRRB 
PRZM,  AGNPS,  HSPF,  WEPP,  and  SWAT.  Several  of  these  models  have 
been  applied  in  both  settings,  but  are  categorized  based  on  their  initial 
development. 


DR3M 


The  Distributed  Routing,  Rainfall,  Runoff  Model  (DR3M)  was  devel¬ 
oped  by  the  U.S.  Geological  Survey  (USGS)  to  simulate  urban  runoff  at 
various  points  in  the  watershed.  The  hydrographs  generated  by  DR3M  are 
used  by  a  companion  model,  DR3M-QUAL,  to  simulate  runoff  quality 
(Alley  and  Smith  1982b). 


Hydrology 

DR3M  (Alley  and  Smith  1982a)  breaks  runoff  into  four  types  of  seg¬ 
ments.  overland,  channel,  reservoir,  and  nodal.  Each  segment  computes 
runoff  differently.  Runoff  from  overland  segments  is  computed  using  the 
kinematic  wave  approximation.  A  daily  soil  moisture  balance  is  main¬ 
tained  using  pan  evaporation.  Infiltration  is  approximated  using  Green- 
Ampt  methodology  to  compute  rainfall  excess.  Two  types  of  overland 
segments  are  considered,  pervious  and  impervious.  A  portion  of  the  rain¬ 
fall  may  be  stored  in  impervious  areas  but  does  not  infiltrate.  The  time 
step  and  distance  used  in  the  kinematic  wave  approximation  are  set  by  the 
user.  Routing  in  the  channel  segments  is  also  done  using  kinematic  wave. 
The  user  must  specify  the  channel  shape  and  size.  Flow  into  a  pipe  may 
be  simulated  by  specifying  a  circular  channel  shape.  Reservoir  routing  can 
be  simulated  by  linear  storage  or  modified  pulse  routing.  Nodal  segments 
are  those  where  three  or  more  of  the  other  segments  meet.  Flow  out  of  the 
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nodal  segment  is  computed  as  the  sum  of  the  segments  entering  the  nodal 
segment.  Nodal  segments  may  also  be  used  to  input  a  hydrograph  or 
remove  runoff. 


Sediment  and  other  pollutants 

The  production  of  sediment  is  computed  in  DR3M-QUAL  (Alley  and 
Smith  1982b)  as  a  function  of  time  between  storms,  build  up/wash  off,  for 
the  impervious  areas.  Sediment  from  the  pervious  areas  is  computed  using 
a  regression  relationship  developed  from  the  USLE.  Other  constituents 
may  be  simulated  as  a  fraction  of  the  sediment.  Soluble  contaminants  are 
not  modeled.  Pollutants  in  the  channels  are  assumed  not  to  mix,  and  the 
transport  is  simulated  using  Lagrangian  methods.  It  is  assumed  reservoirs 
are  plug  flow  and  do  not  mix.  Settling  is  allowed  in  reservoirs  following 
Stokes’  Law. 


Model  use 

DR3M-QUAL  is  not  presently  available  in  PC  form.1  DR3M  is  PC 
compatible,  and  the  source  code  for  the  quality  portion  of  the  model  is 
accessible.  DR3M  has  received  extensive  use  and  review  by  USGS. 
DR3M-QUAL  has  been  used  in-house  by  USGS  to  a  lesser  extent. 


SWMM 


The  Storm  Water  Management  Model  (SWMM)  was  developed  in 
1969-71  for  the  Environmental  Protection  Agency  (EPA)  by  Metcalf  and 
Eddy,  Inc.,  in  conjunction  with  University  of  Florida  and  Water  Resources 
Engineers,  Inc.  SWMM  has  been  updated  several  time  over  the  years  with 
version  4.30  being  released  in  May  1994  (Center  for  Exposure  Assess¬ 
ment  Modeling  (CEAM)  1994).  SWMM  is  one  of  the  most  widely  used 
urban  water  quality  models  (Wurbs  1995). 


Hydrology 

Runoff  and  routing  in  SWMM  are  more  hydraulic  than  hydrologic 
(Viessman  et  al.  1977).  Runoff  from  small  subcatchments  (single  parking 
lots,  city  lots,  etc.)  is  first  routed  overland  to  gutters  or  storm  drains  using 
Manning’s  equation  and  continuity.  Infiltration  is  computed  using  Hor¬ 
ton’s  equation  or  Green- Ampt  methodology.  Depression  storage  fills 
before  overland  flow  begins.  Evapotranspiration  is  neglected.  Once  the 
overland  flow  enters  a  gutter  or  drain,  the  flow  is  combined  with  upstream 


1  Personal  Communication,  1995,  Kate  Flynn,  USGS,  Reston,  VA. 
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flow  and  routed  to  the  final  outflow.  The  Extran  Block  in  SWMM  solves 
the  complete  dynamic  flow  routing  equations  (St.  Venant’s  equation)  to 
accurately  simulate  backwater,  looped  connections,  surcharging,  and  pres¬ 
sure  flow  (CEAM  1994).  SWMM  may  easily  be  calibrated  using  measured 
hydrographs  at  any  point  in  the  flow  system. 


Sediment  and  other  pollutants 

Several  methods  are  available  to  simulate  water  quality.  Sediment  pro¬ 
duction  from  pervious  areas  is  accomplished  using  USLE.  In  impervious 
areas,  build-up/wash-off  methods,  rating  curves,  or  constant  concentration 
can  be  used.  SWMM  also  simulates  deposition  and  scour  in  sewers.  First- 
order  decay  of  constituents  may  be  simulated  (CEAM  1994). 


Model  use 


SWMM  is  widely  used  for  urban  storm  water  planning  and  has  been 
thoroughly  validated.  For  this  purpose,  SWMM  is  an  excellent  model; 
however,  it  cannot  be  applied  readily  to  rural  watersheds  (Viessman  et  al. 
1977).  SWMM  has  been  continually  updated  and  is  supported  by  EPA’s 
CEAM.  A  commercial  PC  version,  XP-SWMM,  is  available  with  better 
data  input  interfaces. 


STORM 


In  the  early  1970’s  the  Hydrologic  Engineering  Center  (HEC)  devel¬ 
oped  the  Storage,  Treatment,  Overflow  Runoff  Model  (STORM)  for  use  in 
urban  areas.  The  output  from  STORM  is  divided  into  statistics  on  runoff 
quality  and  quantity  and  pollutographs  for  individual  events.  The  main 
use  of  STORM  has  been  for  sizing  control  structures  (HEC  1977).  Use  of 
STORM  has  decreased  in  recent  years,  and  HEC  no  longer  supports  the 
model;  however,  several  private  firms  have  updated  STORM  and  support 
it  as  ProSTORM.1 


Hydrology 

STORM  (HEC  1977)  has  three  methods  to  compute  surface  runoff,  the 
SCS  curve  number  method,  a  runoff  coefficient  method,  or  a  combination 
method  where  curve  numbers  are  used  on  pervious  areas  and  coefficients 
are  used  on  impervious  areas.  Evaporation  and  rainfall  are  user  inputs. 
STORM  computes  soil  moisture  conditions  from  evaporation  data  to 
select  the  proper  curve  number.  Dry  weather  flows  are  provided  by  the 


1  Personal  Communication,  1995,  Arlen  Feldman,  HEC,  Davis,  CA. 
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user  to  simulate  combined  sewers.  Flow  is  routed  using  a  unit  hydrograph 
method.  Storage,  treatment,  and  overflow  are  computed  hourly  during  rain¬ 
fall  events,  but  not  during  dry  periods. 


Sediment  and  other  pollutants 

Sediment  production  from  pervious  areas  is  computed  using  the  USLE. 
Sediment  production  from  other  areas  is  based  on  daily  build-up/wash-off 
relationships.  Nitrogen,  phosphorus,  biochemical  oxygen  demand  (BOD), 
and  coliforms  may  also  be  simulated,  but  they  are  only  represented  by 
loading  functions  and  do  not  degrade  or  cycle  (HEC  1977). 


Model  use 

STORM  was  applied  in  several  runoff  studies  in  the  late  1970’s  and 
early  1980’s  (Donigian  and  Huber  1990).  Abbott  (1978)  compared 
STORM  to  other  continuous  and  event  models  as  well  as  actual  data  and 
achieved  good  results.  HEC  no  longer  supports  STORM,  and  a  PC  version 
is  currently  available.  ProSTORM  is  available  for  PC  use  from  private 
consulting  firms.1 


CREAMS/GLEAMS 


The  Chemical,  Runoff,  and  Erosion  from  Agricultural  Management  Sys¬ 
tems  (CREAMS)  model  was  the  first  in  a  series  of  models  developed  by 
the  U.S.  Department  of  Agriculture  (USD A).  The  original  model  focused 
on  surface  water  quality,  but  groundwater  quality  was  added  to  form  the 
Groundwater  Loading  Effects  of  Agricultural  Management  Systems 
(GLEAMS)  model  (Knisel  1993).  GLEAMS  has  replaced  CREAMS  for 
the  most  part,  but  CREAMS  is  the  recognized  predecessor  for  several 
other  models.  CREAMS  is  a  field-scale  model  used  to  look  at  edge  of 
field  loading  sediment,  nutrients,  and  pesticides  with  various  management 
practices.  With  the  additions  in  GLEAMS,  leaching  below  the  root  zone 
may  also  be  examined  (Knisel  1993).  CREAMS/GLEAMS  is  one  of  the 
most  detailed  field-scale  models  currently  in  use.  USDA’s  Agriculture 
Research  Service  (USDA-ARS)  has  recently  developed  the  Soil  and  Water 
Assessment  Tool  (SWAT)  by  extending  GLEAMS  to  basin  scale. 


1  Personal  Communication,  1995,  Arlen  Feldman,  HEC,  Davis,  CA, 
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Hydrology 

The  hydrology  component  of  the  model  is  based  on  a  daily  water  bal¬ 
ance.  Surface  runoff  is  based  on  the  SCS  curve  number  method  when 
daily  rainfall  data  are  used.  If  breakpoint  rainfall  data  are  used,  runoff  is 
computed  using  Green-Ampt  infiltration  and  kinematic  wave  routing.  The 
peak  runoff  rate  is  estimated  by  regression  relationships  (Smith  and  Wil¬ 
liams  1980).  Percolation  is  computed  using  a  storage-routing  technique  to 
predict  flow  through  each  soil  layer.  Lateral  subsurface  flow  is  computed 
for  each  soil  layer  using  a  kinematic  storage  model  starting  at  the  top 
layer  and  progressing  downward  (Knisel  1993).  Soil  evaporation  may  be 
computed  using  the  Priestly-Taylor  method  or  Penman-Monteith.  Plant 
evapotranspiration  is  based  on  the  Ritchie  method  and  requires  the  leaf 
area  index  depending  on  the  stage  of  crop  growth.  Snowmelt  and  transmis¬ 
sion  losses  are  also  taken  into  account  (Knisel  1993).  The  CREAMS/ 
GLEAMS  models  have  no  channel  routing  since  they  are  field-scale 
models. 


Sediment  and  other  pollutants 

Sediment  yield  in  CREAMS  is  computed  using  detachment,  deposition, 
and  transport  equations.  Detachment  is  divided  into  two  categories,  rill 
and  interrill.  Detachment  in  the  interrill  areas  is  based  on  rainfall  energy, 
and  rill  erosion  is  based  on  flow  velocities.  The  interrill  erosion  is  com- ' 
puted  first;  then  the  transport  capacity  is  computed.  If  the  transport  capac¬ 
ity  of  the  overland  flow  exceeds  the  interrill  sediment  load,  rill  erosion  is 
computed.  The  sediment  load  is  the  lesser  of  transport  capacity  and  the 
combined  rill-interrill  sediment  load.  If  the  interrill  sediment  load 
exceeds  transport  capacity,  deposition  occurs.  Five  particle  sizes  are  con¬ 
sidered  for  sediment  transport  and  detachment/deposition:  sand,  silt,  clay, 
small  aggregates,  and  large  aggregates.  Deposition  in  small  impound¬ 
ments  and  transport  and  deposition  in  terraces  or  small  channels  is  also 
considered  (Foster  et  al.  1980). 

CREAMS/GLEAMS  has  a  well-documented  submodel  for  dealing  with 
pesticides  that  can  handle  up  to  20  pesticides  at  once  (Knisel  1993).  The 
amount  of  pesticide  lost  to  runoff  is  computed  as  a  soluble  fraction  and  an 
adsorbed  fraction  of  the  available  pesticide  that  has  not  leached,  volatil¬ 
ized,  decayed,  or  otherwise  been  lost.  Different  methods  of  pesticide  appli¬ 
cation,  foliar  or  soil  incorporated,  are  considered,  and  the  effects  on  the 
initial  concentrations  of  pesticides  are  computed.  Using  GLEAMS,  the 
amount  of  pesticide  that  leaches  to  the  groundwater  may  also  be  deter¬ 
mined  (Knisel  1993).  Since  the  model  only  looks  at  edge  of  field  concen¬ 
tration,  no  degradation  is  considered  in  the  channels. 

Nitrogen  losses  are  considered  for  both  nitrate  and  organic  nitrogen. 

For  nitrate,  the  daily  soil  nitrate  content  is  computed  for  each  soil  layer. 

The  amount  of  nitrate  that  is  leached  from  each  layer  is  computed  and 
will  percolate  to  the  groundwater  or  may  reappear  with  return  flow. 
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Nitrate  in  the  runoff  is  computed  using  the  nitrate  in  the  top  soil  layer 
only.  Organic  nitrogen  is  considered  with  the  sediment  with  a  loading 
function.  Crop  uptake  of  nitrogen  is  also  considered.  The  different  proc¬ 
esses  that  affect  nitrogen  are  considered,  mineralization,  immobilization, 
denitrification,  volatilization,  and  fixation  (Knisel  1993).  Phosphorus  is 
mainly  transported  with  sediment  such  as  organic  nitrogen.  Soluble  phos¬ 
phorus  transported  by  runoff  is  computed  from  the  level  of  phosphorus  in 
the  top  soil  layer.  Application  of  fertilizers  and  crop  growth  and  residue 
are  taken  into  account  when  computing  the  nutrient  levels  in  each  soil 
layer  (Knisel  1993). 


Model  use 

The  CREAMS/GLEAMS  models  have  been  widely  used  in  many  areas 
(Knisel  1993).  Long-term  studies  on  well-instrumented  fields  have  shown 
that  the  methodologies  used  in  CREAMS/GLEAMS  work  well.  The  mod¬ 
els  are  readily  available  for  PC  use. 


EPIC 


Concerns  over  the  possible  decrease  in  agricultural  productivity  due  to 
soil  loss  prompted  the  creation  of  the  Erosion/Productivity  Impact  Calcu¬ 
lator  (EPIC)  model.  EPIC  was  developed  to  determine  the  long-term  pro¬ 
ductivity  of  soil  using  various  management  strategies  (Williams  1994  and 
Dumesnil  1993).  Simulations  using  EPIC  are  usually  done  on  small  areas 
over  a  long  period  of  time.  Since  the  model’s  main  purpose  is  to  compute 
productivity,  a  greater  portion  of  the  model  deals  with  crop  growth  com¬ 
pared  with  other  models,  but  water  yield,  sediment  production,  and  nutri¬ 
ent  removal  are  also  considered  (Williams  1994). 


Hydrology 

The  hydrology  component  of  the  model  is  based  on  a  daily  water  bal¬ 
ance.  Surface  runoff  is  based  on  the  SCS  curve  number  method,  and  the 
peak  runoff  rate  is  estimated  by  the  rational  formula.  Percolation  is  com¬ 
puted  using  a  storage-routing  technique  to  predict  flow  through  each  soil 
layer.  Once  water  has  percolated  past  the  root  zone,  it  is  considered 
groundwater  and  lost  from  the  model.  Lateral  subsurface  flow  is  com¬ 
puted  for  each  soil  layer  using  a  kinematic  storage  model  starting  at  the 
top  layer  and  progressing  downward.  Soil  evaporation  may  be  computed 
using  either  the  Penman-Monteith  method  or  the  Priestly-Taylor  method. 
Plant  evapotranspiration  is  based  on  the  Ritchie  method  and  requires  the 
leaf  area  index  depending  on  the  stage  of  crop  growth.  Snowmelt  and 
transmission  losses  are  also  taken  into  account.  Since  EPIC  is  a  small- 
scale  model,  channel  routing  is  not  considered.  Irrigation  of  crops  may  be 
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considered  with  the  daily  water  balance  (Williams  1994  and  Sharpley  and 
Williams  1990). 


Sediment  and  other  pollutants 

Sediment  yield  may  be  computed  in  one  of  three  ways  specified  by  the 
user,  USLE,  Modified  USLE  (MUSLE),  or  the  Onstad-Foster  modification 
of  the  USLE.  Erosion  from  furrow  irrigation  is  always  computed  using 
MUSLE.  The  variation  between  these  models  is  the  energy  factor  used  to 
drive  erosion;  USLE  uses  rainfall  only,  MUSLE  uses  runoff  only,  and 
Onstad-Foster  uses  a  combination  of  rainfall  and  runoff.  As  sediment 
leaves  the  field,  the  surface  layer  of  soil  is  reduced  in  thickness  until  it  no 
longer  remains.  Erosion  continues  to  move  through  the  soil  profile,  allow¬ 
ing  for  soil  weathering  (Williams  1994  and  Sharpley  and  Williams,  1990). 
Wind  erosion  is  computed,  but  washoff  does  not  affect  runoff  sediment 
loads  (Sharpley  and  Williams  1990). 

Losses  of  nitrate  are  considered  from  the  top  layer  of  soil  only.  Nitrate 
that  is  adsorbed  and  in  solution  may  leave  with  the  runoff.  Nitrate  in  solu¬ 
tion  may  also  leave  the  field  though  percolation  or  lateral  subsurface  flow. 

Loading  functions  for  organic  nitrogen  are  also  provided.  Nitrogen 
may  be  moved  into  the  top  soil  layer  as  water  in  that  layer  evaporates  and 
water  from  lower  layers  enters.  The  nitrogen  cycle  is  simulated  by  the 
processes  of  denitrification,  mineralization,  immobilization,  and  fixation. 
The  concentration  of  nitrogen  in  rainfall  may  also  be  included.  Phospho¬ 
rus  is  transported  in  solution  or  with  sediment.  The  phosphorus  cycle 
includes  mineralization,  immobilization,  and  mineral  cycling.  EPIC  does 
not  track  pesticide  losses  (Williams  1994,  and  Sharpley  and  Williams 
1990). 


Other  features 

EPIC  is  able  to  generate  long-term  sediment  predictions  because  of  its 
built-in  weather  generator  and  weather  database  (Sharpley  and  Williams 
1990).  Rainfall  timing  is  based  on  the  probability  of  a  wet  day  following  a 
dry  day  and  a  wet  day  following  a  wet  day.  Daily  rainfall  amounts  are 
predicted  using  a  skewed  normal  distribution  with  the  mean  rainfall  and 
standard  deviation  for  the  month.  Air  temperature  and  solar  radiation  are 
correlated  with  rainfall.  Wind  speed  is  predicted  from  a  two-parameter 
gamma  distribution,  and  direction  is  based  on  a  cumulative  probability  dis¬ 
tribution.  Daily  relative  humidity  is  derived  using  the  monthly  average 
and  a  triangular  distribution  and  adjusting  for  a  wet  or  dry  day.  Soil  tem¬ 
perature  and  soil  pH  are  also  computed  and  applied  to  the  nutrient  cycle. 


Chapter  2  Model  Review 


EPIC  includes  a  crop  growth  simulator  to  account  for  various  stages  of 
crop  growth  and  their  impacts  such  as  nutrient  uptake,  transpiration,  and 
soil  residue.  An  economic  section  is  included  to  further  assess  soil  produc¬ 
tivity  (Williams  1994,  and  Sharpley  and  Williams  1990). 


Model  use 

EPIC  has  been  used  to  simulate  many  different  things,  soil  productiv¬ 
ity,  crop  growth,  and  soil  degradation  (Dumesnil  1993).  Sharpley  and 
Williams  (1990)  report  on  the  validation  of  the  model  components  and 
application  of  the  complete  model  to  differing  soils  and  climates.  Sensitiv¬ 
ity  analysis  was  also  reported.  EPIC  is  available  for  PC  use  and  is  well 
documented  and  supported  by  USDA-ARS. 


SWRRB 

The  SWRRB  (Simulator  for  Water  Resources  in  Rural  Basins)  model 
was  developed  by  the  USDA-ARS  to  determine  the  effects  of  management 
practices  on  water  and  sediment  yields  in  small  ungaged  rural  basins 
(Arnold  and  Williams  1994).  The  model  was  developed  by  the  USDA 
extending  the  CREAMS/GLEAMS  methodology  from  field  scale  to  basin 
scale.  The  major  processes  included  in  the  model  are  surface  runoff,  perco¬ 
lation,  return  flow,  transmission  losses,  pond  and  reservoir  storage,  sedi¬ 
mentation,  pesticide  fate,  nutrient  cycling,  and  crop  growth.  Several 
subbasins  with  varying  soil,  crop,  and  weather  conditions  can  be  simu¬ 
lated  as  part  of  a  single  watershed.  The  model  was  designed  to  look  at 
long-term  effects  of  management  strategies,  but  single  runoff  events  can 
be  simulated. 


Hydrology 

The  hydrology  component  of  the  model  is  based  on  a  daily  water  bal¬ 
ance.  Surface  runoff  is  based  on  the  SCS  curve  number  method,  and  the 
peak  runoff  rate  can  be  estimated  by  either  the  rational  formula  or  SCS 
TR-55  method.  Percolation  is  computed  using  a  storage-routing  technique 
with  a  crack  flow  model  to  predict  flow  through  each  soil  layer.  Once 
water  has  percolated  past  the  root  zone,  it  is  considered  groundwater,  but 
may  appear  as  return  flow  in  downstream  basins.  Lateral  subsurface  flow 
is  computed  for  each  soil  layer  using  a  kinematic  storage  model  starting  at 
the  top  layer  and  progressing  downward.  Soil  evaporation  may  be  com¬ 
puted  using  either  the  Hargreaves-Samani  method  or  the  Priestly-Taylor 
method.  Plant  evapotranspiration  is  based  on  the  Ritchie  method  and 
requires  the  leaf  area  index  depending  on  the  stage  of  crop  growth.  Snow¬ 
melt  and  transmission  losses  are  also  taken  into  account  (Arnold  and 
Williams  1994).  Since  SWRRB  is  a  long-term  water  and  sediment  yield 
predictor,  flood  routing  is  not  performed.  Daily  basin  outflow  is  computed 
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by  summing  the  subbasin  outflow.  Farm  ponds  and  reservoirs  can  be  simu¬ 
lated  with  SWRRB.  Irrigation  of  crops  may  be  considered  with  the  daily 
water  balance  (Arnold  and  Williams  1994). 


Sediment  and  other  pollutants 

Sediment  yield  is  computed  for  each  subbasin  using  MUSLE.  Sediment 
in  channel  routing  consists  of  two  processes  operating  simultaneously, 
deposition  and  degradation.  Stokes’  Law  is  applied  to  each  sediment  size 
for  each  channel  reach  to  determine  deposition.  Channel  degradation  is  a 
function  of  the  stream  power.  Sediment  is  allowed  to  be  deposited  and 
reentrained  in  the  channel  before  reaching  the  basin  outlet.  A  sediment 
balance  is  maintained  for  ponds  and  reservoirs  as  well  (Arnold  and 
Williams  1994). 

SWRRB  uses  the  same  methodology  as  GLEAMS  (Knisel  1993)  to 
deal  with  pesticides.  First  the  amount  of  pesticide  reaching  the  ground  is 
computed  dependent  on  the  stage  of  crop  growth.  The  amount  of  pesticide 
that  leaches  to  the  groundwater  is  computed.  The  amount  of  pesticide  lost 
to  runoff  is  computed  as  a  soluble  fraction  and  an  adsorbed  fraction  of  the 
available  pesticide  that  has  not  leached,  volatilized,  decayed,  or  otherwise 
been  lost.  Pesticides  that  are  adsorbed  to  sediment  are  allowed  to  be 
deposited  in  the  channel.  Soluble  pesticides  are  considered  conservative 
in  the  channel  and  do  not  degrade  (Arnold  and  Williams  1994). 

Nitrogen  losses  are  considered  for  both  nitrate  and  organic  nitrogen. 

For  nitrate,  the  daily  soil  nitrate  content  is  computed  for  each  soil  layer. 
The  amount  of  nitrate  that  is  leached  from  each  layer  is  computed  and 
will  percolate  to  the  groundwater  or  may  reappear  with  return  flow. 

Nitrate  in  the  runoff  is  computed  using  the  nitrate  in  the  top  soil  layer 
only.  Organic  nitrogen  is  considered  with  the  sediment  with  a  loading 
function.  Crop  uptake  of  nitrogen  is  also  considered.  Phosphorus  is 
mainly  transported  with  sediment  such  as  organic  nitrogen.  Soluble  phos¬ 
phorus  transported  by  runoff  is  computed  from  the  level  of  phosphorus  in 
the  top  soil  layer,  but  leaching  is  not  considered.  Soluble  nutrients  are  con¬ 
sidered  conservative  in  the  basin  channels,  but  adsorbed  nutrients  may  be 
deposited  or  entrained  in  the  channel.  Application  of  fertilizers  and  crop 
growth  and  residue  are  taken  into  account  when  computing  the  nutrient 
levels  in  each  soil  layer  (Arnold  and  Williams  1994). 


Other  features 

SWRRB  allows  the  user  to  input  the  required  weather  data  (daily  rain¬ 
fall  and  temperature),  or  a  stochastic  weather  generator  can  simulate  the 
weather.  Solar  radiation  is  always  simulated.  The  probability  of  a  rainfall 
event  is  simulated  using  a  first-order  Markov  chain.  When  rainfall  does 
occur,  the  amount  is  simulated  from  a  skewed  normal  distribution.  Tem¬ 
perature  and  solar  radiation  are  simulated  using  Richardson’s  model 
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which  correlates  temperature  and  radiation  with  precipitation.  Weather 
may  be  simulated  for  the  entire  basin  or  each  subbasin  individually.  The 
parameters  needed  for  the  weather  generator  are  included  with  the  model 
for  locations  across  the  United  States  including  Puerto  Rico  (Arnold  and 
Williams  1994). 

Tillage  practices,  crop  residue,  and  crop  growth  are  simulated  as  well 
as  how  they  affect  soil  properties.  Data  for  numerous  soils  that  may  be 
used  are  included  with  the  model  or  may  be  input  by  the  user  for  each  sub¬ 
basin.  A  utility  program  is  included  to  assist  with  data  entry,  but  the  input 
is  still  in  tabular  form.  The  output  is  in  tabular  form.  SWRRB  has  been 
combined  successfully  with  a  lake  water  quality  model  to  simulate  the  end 
effects  of  management  decisions  (Arnold  and  Williams  1994). 


Model  use 

SWRRB  has  been  used  and  validated  in  a  variety  of  settings.  The 
model  has  been  used  to  predict  urbanization,  coastal  pollution,  pesticide 
assessment,  pond  irrigation,  and  water  rights  (Arnold  and  Williams  1994). 
Support  and  documentation  for  the  PC  version  of  SWRRB  are  available 
from  USDA-ARS. 


WEPP 

The  objective  of  the  Water  Erosion  Prediction  Project  (WEPP)  is  to 
develop  a  new-generation,  process-oriented  model  to  predict  soil  erosion. 
The  WEPP  model  represents  the  start  of  a  new  generation  of  USDA  mod¬ 
els.  It  was  developed  for  continuous  simulations  of  small  watersheds,  less 
than  260  ha,  or  individual  hillslopes.  WEPP  allows  the  user  to  look  at  spa¬ 
tial  and  temporal  variations  in  sediment  production.  While  this  newly 
developed  model  uses  some  of  the  most  up-to-date  methodology,  since 
some  components  are  missing  it  is  not  a  complete  water  quality  model 
(Flanagan  et  al.  1995). 


Hydrology 

The  hydrology  of  WEPP  is  more  complex  and  physically  based  than 
SWRRB.  The  hydrology  component  of  the  model  is  based  on  a  daily 
water  balance.  Surface  runoff  is  computed  as  the  rainfall  excess,  rainfall 
less  infiltration,  rather  than  using  an  empirical  method.  Infiltration  is  com¬ 
puted  using  the  Green-Ampt  equation  as  modified  by  Mein  and  Larson. 
Peak  runoff  and  runoff  rate  are  computed  by  routing  the  rainfall  excess 
downslope  using  the  kinematic  wave  approximation.  Percolation  is  com¬ 
puted  using  methods  similar  to  SWRRB  as  well  as  simulating  flow  to 
drain  tiles  and  ditches.  WEPP  computes  soil  evaporation  using  either  the 
Penman  method  or  the  Priestly-Taylor  method.  Plant  evapotranspiration  is 
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based  on  the  Ritchie  method  and  requires  the  leaf  area  index  depending  on 
the  stage  of  crop  growth  (Stone  et  al.  1995).  Snowmelt  and  transmission 
losses  are  also  taken  into  account  (Savabi  et  al.  1995).  In  the  watershed 
mode,  WEPP  routes  hillslope  runoff  through  channels  to  compute  the 
basin  outflow.  Farm  ponds  and  reservoirs  can  be  simulated  as  with 
SWRRB  (Ascough  et  al.  1995).  Irrigation  of  crop  may  be  considered  with 
the  daily  water  balance. 


Sediment  and  other  pollutants 

WEPP  computes  sediment  yield  in  a  method  similar  to  CREAMS/ 
GLEAMS.  Sediment  production  is  divided  into  two  parts,  rill  and  interrill. 
Detachment  in  the  interrill  areas  is  a  function  of  rainfall  intensity,  and  the 
detachment  process  is  dominated  by  raindrop  impact.  All  interrill  sedi¬ 
ment  is  carried  to  the  rills  or  small  channels.  The  transport  capacity  of  the 
rill  is  then  computed,  and  interrill  sediment  is  either  deposited  or  trans¬ 
ported.  If  the  rill  has  an  excess  transport  capacity,  detachment  of  sediment 
in  the  rill  is  assumed  to  occur.  Various  soil  particle  sizes  are  considered 
(Foster  et  al.  1995).  Unlike  CREAMS/GLEAMS,  transport  capacity  and 
shear  stress  is  based  on  rill  hydraulics  rather  than  sheet  flow.  Sediment  is 
allowed  to  be  deposited  and  detached  in  the  channel  when  using  the  water¬ 
shed  version  of  WEPP  (Ascough  et  al.  1995).  Unlike  the  other  model, 
WEPP  does  not  rely  on  the  USLE  to  incorporate  management  factors,  but 
rather  focuses  on  the  physical  processes. 


Other  features 

WEPP  uses  the  weather  generator  developed  for  EPIC  and  SWRRB 
with  some  additions  and  refinements  (Nicks,  Lane,  and  Gander  1995).  The 
probability  of  a  rainfall  event  is  simulated  using  a  first-order  Markov 
chain.  When  rainfall  does  occur,  the  amount  is  simulated  from  a  skewed 
normal  distribution.  Because  WEPP  does  not  use  the  SCS  method  to  com¬ 
pute  runoff,  storm  duration,  peak  intensity,  and  break  point  rainfall  are 
computed.  Temperature  and  solar  radiation  are  simulated  using  Richard¬ 
son  s  model  which  correlates  temperature  and  radiation  with  precipitation. 
Dew  point,  wind  speed,  and  direction  are  additions  for  use  by  the  evapo- 
transpiration  computations.  Weather  may  be  simulated  for  the  entire  basin 
or  each  hillslope  individually.  The  weather  files  included  with  WEPP  are 
not  as  extensive  as  SWRRB  at  this  time.  Linking  WEPP  with  GIS  and  the 
1995°naI  ^^mat'C  ^ata  Center  is  currently  being  explored  (Flanagan  et  al. 

Tillage  practices,  crop  residue,  and  crop  growth  are  simulated  as  well 
as  how  they  affect  soil  properties.  Data  for  several  soils  that  may  be  used 
are  included  with  the  model  or  may  be  input  by  the  user  for  each  sub¬ 
basin.  The  user  can  use  several  previously  developed  files  to  create  a 
watershed  by  combining  soils,  weather,  conservation  practices,  slope,  and 
channel  files. 
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Model  use 


Many  of  the  components  of  WEPP  have  been  validated  separately,  but 
no  widespread  use  of  WEPP  has  been  reported.  The  model  is  still  under 
development  and  has  yet  to  be  fully  verified  and  validated.  Components 
of  WEPP  may  be  used  in  other  models.  PC  versions  of  the  model,  as  well 
as  documentation,  support,  and  development,  are  available  from  the 
USDA-ARS  National  Soil  Erosion  Laboratory. 


SWAT 

The  Soil  and  Water  Assessment  Tool  (SWAT)  model  is  the  latest  model 
developed  by  the  USDA-ARS.  Where  SWRRB  was  an  extension  of 
CREAMS/GLEAMS  from  field  scale  to  small  watershed  scale,  SWAT  is 
the  extension  of  SWRRB  to  basin  scale  (Arnold  et  al.  1995).  SWAT  was 
developed  by  including  ROTO  (Routing  Output  to  Outlets)  with  SWRBB. 
Complete  flow  and  contaminant  routing  is  accomplished  as  well  as  the 
inclusion  of  GIS  interfaces,  a  weather  generator,  lake  water  quality,  and 
water  management  options. 


Hydrology 

Runoff  is  estimated  using  daily  rainfall  data  using  the  SCS  curve 
number  method  with  peak  runoff  computed  using  a  modification  of  the 
rational  method.  Time  to  concentration  for  both  channels  and  overland 
flow  is  estimated  using  Manning’s  equation.  Percolation  and  lateral  sub¬ 
surface  flow  are  computed  together  and  are  based  on  storage  routing  and 
kinematic  routing.  Groundwater  flow  is  simulated  based  on  a  shallow  aqui¬ 
fer  model  with  evaporation,  pumping,  seepage,  and  discharge  being 
included.  The  daily  water  balance  includes  evapotranspiration  based  on 
Hargreaves-Samani,  Priestly-Taylor,  or  Penman-Monteith  equations.  Plant 
transpiration  is  based  on  the  Ritchie  method.  Transmission  losses  in  chan¬ 
nels  are  considered  along  with  ponds  and  other  small  structures  that  will 
reduce  runoff  quantity  (Arnold  et  al.  1995).  Water  may  be  transferred 
from  one  portion  of  the  watershed  to  another  to  simulate  pumping  and 
irrigation. 


Sediment  and  other  pollutants 

Sediment  yield  is  based  on  the  MUSLE.  Nitrogen,  phosphorus,  and  pes¬ 
ticides  are  simulated  using  the  CREAMS/GLEAMS  (Knisel  1993)  models 
including  biochemical  processes  and  groundwater  loading  as  described 
earlier  in  this  report.  During  runoff  events,  chemical  degradation  is  not 
simulated  in  the  channels,  but  a  lake  water  quality  model  has  been  added 
to  simulate  lake  processes.  Sediment-bound  chemicals  may  settle  and  be 
detached  in  channels.  Crop  growth  simulators  are  included  to  account  for 
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changing  ground  cover  and  plant  uptake  of  nutrients.  Point  sources  may 
also  be  included  in  the  simulations  (Arnold  et  al.  1995). 


Other  features 

SWAT  includes  a  weather  generator  similar  to  the  one  used  in  EPIC, 
SWRRB,  and  WEPP  for  long-term  simulations  and  numerous  weather  sta¬ 
tion  locations  to  reduce  user  input  (Arnold  et  al.  1995).  Weather  may  be 
simulated  for  the  entire  basin  or  for  each  subwatershed  individually  as  is 
done  in  SWRRB  (Arnold  and  Williams  1994).  SWAT  also  includes  an 
interface  with  GRASS  (Geographic  Resources  Analysis  Support  System) 
GIS  (Arnold  et  al.  1995).  This  allows  the  user  to  easily  input  data  from  a 
large  heterogeneous  basin.  One  feature  in  SWAT  not  seen  in  any  other 
model  is  the  mixing  of  watershed  delineation  techniques.  A  grid  system, 
as  seen  in  AGNPS  (Agricultural  Non-Point  Source  pollution  model),  may 
be  used  in  areas  of  special  interest  while  the  rest  of  the  basin  may  be  simu¬ 
lated  as  more  homogeneous  subwatersheds.  This  allows  detailed  simula¬ 
tion  of  areas  of  interest  without  a  large  increase  in  input  data  (Arnold 
et  al.  1995). 


Model  use 


SWRRB  and  ROTO  have  both  been  extensively  tested.  SWAT  has  been 
applied  to  several  gaged  watersheds  and  rendered  good  results.  PC  ver¬ 
sions  of  SWAT  along  with  support  and  documentation  are  available  from 
USDA-ARS. 


HSPF 


The  Hydrologic  Simulation  Program — FORTRAN  (HSPF)  model  may 
be  one  of  the  most  complete  models  to  deal  with  urban  and  nonurban  run¬ 
off  (Donigian  and  Huber  1990).  The  hydrology  of  HSPF  is  based  on  the 
Stanford  Watershed  Model,  developed  in  the  early  1960’s.  Two  erosion 
components  deal  with  pervious  and  impervious  areas  separately  to  avoid 
forcing  one  type  of  model  to  operate  in  an  area  it  was  not  designed  to 
simulate.  HSPF  contains  many  submodels  to  simulate  many  of  the  physi¬ 
cal,  chemical,  and  biological  processes  in  a  watershed  (Bicknell  et  al. 
1993). 


Hydrology 

HSPF  is  based  on  the  Stanford  Watershed  Model.  HSPF  uses  basic  rela¬ 
tionships  to  represent  the  physical  processes.  The  model,  like  most,  is 
based  on  the  water  balance  of  the  soil.  Rainfall  rate  is  a  user  input.  Rain¬ 
fall  is  reduced  by  an  interception  rate  based  on  the  type  of  ground  cover. 
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Evapotranspiration  is  based  on  adjusted  pan  evaporation  data.  The  loca¬ 
tion  that  the  evaporation  comes  from  is  based  on  the  available  water  in 
certain  areas.  Infiltration  is  considered  to  vary  linearly  with  time  and  have 
a  maximum  infiltration  capacity.  Interflow  storage  is  also  considered.  In 
the  impervious  areas,  infiltration,  interflow,  and  soil  evaporation  are 
assumed  to  not  occur.  Overland  flow  is  based  on  Darcy’s  law  modified  to 
adjust  for  the  amount  of  surface  storage  during  the  time  interval.  Channel 
routing  and  reservoir  storage  are  accomplished  by  using  the  time-area 
method  and  summing  the  inputs  (Bicknell  et  al.  1993). 


Sediment  and  other  pollutants 

Sediment  production  is  based  on  detachment  by  rainfall  and  transport 
by  the  overland  flow  in  a  pervious  area.  Scouring  of  the  soil  matrix  is  also 
considered.  Three  sizes  of  sediment  are  considered;  sand,  silt,  and  clay 
are  considered  in  streams,  but  not  in  erosion.  Impervious  sediment  load¬ 
ing  is  determined  by  build-up/wash-off  relationships.  Sediment  loading  in 
channels  and  reservoirs  is  based  on  the  sum  of  the  incoming  loads,  and 
settling  is  allowed.  Scour  and  deposition  are  based  on  shear  stress  calcula¬ 
tions.  Channels  and  reservoirs  are  considered  well  mixed  (Bicknell  et  al. 
1993). 

Nitrogen,  phosphorus,  and  three  pesticides  may  be  simulated.  Constitu¬ 
ents  may  move  with  sediment,  on  the  surface  in  solution  or  in  the  ground- 
water  in  solution.  Sorption  may  be  simulated  using  first-order  kinetics  or 
Freundlich  isotherms.  HSPF  also  accounts  for  up  to  two  degradation  prod¬ 
ucts  of  pesticides.  Nutrient  cycles  and  instream  kinetics  are  simulated. 
Decay  and  instream  kinetics  include  hydrolysis,  oxidation,  photolysis, 
volatilization,  biodegradation,  and  generalized  first-order  decay  (Bicknell 
et  al.  1993).  Other  constituents  that  may  be  simulated  include  water  tem¬ 
perature,  dissolved  oxygen,  BOD,  and  algal  populations  (benthic,  phyto¬ 
plankton,  and  zooplankton). 


Model  use 

HSPF  has  been  widely  used  in  both  urban  and  nonurban  areas  and  is 
well  validated  in  these  areas  (Donigian  and  Huber  1990).  The  EPA  sup¬ 
ports  HSPF,  and  extensive  documentation  is  available.  Short  courses  on 
the  use  of  HSPF  are  regularly  provided  as  well.  PC  versions  of  the  model 
are  available  and  upgrades  continue.  The  model  is  supported  by  both  EPA 
and  USGS. 
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AGNPS 


Another  model  that  was  developed  for  agricultural  watersheds  is  the 
Agricultural  Non-Point  Source  (AGNPS)  pollution  model.  Unlike  most 
other  models,  AGNPS  divides  the  watershed  into  cells  based  on  a  grid 
pattern.  Each  cell  is  modeled  simultaneously  with  each  other  cell.  The 
ability  to  set  various  parameters  individually  for  each  cell  means  better 
spatial  description  but  more  data  input.  The  model  was  originally  devel¬ 
oped  to  determine  the  impact  of  feedlots  on  runoff  quality,  but  it  may  oper¬ 
ate  with  multiple  point  sources  or  as  a  nonpoint  source  model.  The  model 
will  accommodate  up  to  28,000  cells  ranging  in  size  from  0.1  to  100  acres 
(0.0405  to  40.5  ha),  but  it  is  most  commonly  used  for  smaller  watersheds 
(Young  et  al.  1994). 


Hydrology 

The  current  version  of  AGNPS  is  an  event  model;  therefore,  many  of 
the  hydrologic  components  normally  found  in  other  models  are  missing. 
Runoff  is  computed  using  the  SCS  curve  number  method.  Peak  runoff 
may  either  be  computed  using  SCS  TR-55  or  the  CREAMS  method.  The 
user  must  route  the  flow  from  cell  to  cell.  Storm  duration,  intensity,  and 
type  must  be  entered  for  the  entire  watershed.  The  shape  of  the  runoff 
hydrograph  may  also  be  specified  by  the  user.  Once  the  runoff  leaves  a 
cell  it  is  considered  to  flow  in  channels  to  the  outlet  of  the  watershed.  The 
user  may  specify  channel  shapes  and  flow  characteristics.  The  user  may 
include  impoundments  in  the  channel  system,  and  cells  may  have  no  flow 
outlet  to  simulate  features  such  as  sink  holes  (Young  et  al.  1994). 


Sediment  and  other  pollutants 

Sediment  production  from  each  cell  is  computed  using  the  MUSLE. 
Differing  slopes,  slope  shapes,  soils,  and  management  practices  can  be 
considered  for  each  cell.  When  the  sediment  reaches  a  channel,  the  trans¬ 
port  capacity  of  the  channel  is  computed  for  differing  sediment  sizes  to 
determine  if  there  will  be  deposition  or  channel  detachment.  The  Einstein 
approach  is  used  to  determine  transport  capacity  (Young  et  al.  1994). 
Nitrogen,  phosphorus,  pesticide,  and  chemical  oxygen  demand  concentra¬ 
tions  my  be  tracked  from  either  point  or  nonpoint  sources.  For  point 
sources  the  flow  rate  from  the  cell  and  the  pollutant  concentrations  must 
be  specified.  In  the  case  of  a  feedlot,  runoff  will  be  computed  and  pollut¬ 
ant  loads  may  either  be  specified  or  computed  by  the  model  based  on  the 
size  and  type  of  feedlot.  Nonpoint  pollution  loads  are  calculated  based  on 
the  runoff  and  sediment  loads.  Nutrients  are  allowed  to  decay  in  impound¬ 
ments  and  channels.  The  nitrogen  concentration  of  the  rainfall  may  be 
specified  by  the  user.  Fertilizer  application  rates,  timing,  and  methods 
may  be  specified.  If  a  pesticide  is  simulated,  its  concentration  in  solution 
and  in  the  sediment  is  kept  in  equilibrium  using  its  partition  coefficient. 
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At  this  time  only  one  pesticide  may  be  simulated  for  each  event  (Young 
et  al.  1994). 


Other  features 

The  main  drawback  to  a  grid  type  model  like  AGNPS  is  the  large 
amount  of  data  needed  to  simulate  a  large  watershed.  To  help  alleviate 
this  drawback,  AGNPS  allows  the  use  to  set  larger  cell  sizes  and  then  sub¬ 
divide  the  cell  up  to  64  times  to  simulate  specific  small  features  without 
increasing  the  data  requirements  for  the  entire  watershed.  Onscreen  help 
is  available  for  data  entry  along  with  a  routine  to  check  the  data.  Work  is 
under  way  to  link  AGNPS  with  ARC-INFO  and  GRASS  GIS  databases  to 
assist  with  data  input  and  output  analyses  (Young  et  al.  1994). 


Model  use 

AGNPS  is  available  in  a  PC  version  and  is  well  documented.  Limited 
support  is  available  through  USDA-ARS  (Young  et  al.  1994).  Several  stud¬ 
ies  have  used  AGNPS  to  determine  pollution  loading  from  events  and 
assist  in  the  reduction  of  runoff  pollution  from  feedlots.  The  current  ver¬ 
sion  of  AGNPS  is  an  event  model,  but  work  is  continuing  on  producing  a 
continuous  version.  Research  on  AGNPS  has  been  ongoing  for  numerous 
years. 


PRZM 


The  Pesticide  Root  Zone  Management  (PRZM)  model  was  developed 
to  examine  the  leaching  of  pesticides  through  the  root  zone.  PRZM  has 
been  well  documented  and  validated  for  the  evaluation  of  pesticides.  The 
model  tracks  water  and  pesticide  movement  in  the  vertical  direction  only 
on  a  daily  time  step  (Mullins  et  al.  1996).  Recent  releases  of  PRZM  now 
include  limited  nitrogen  processes. 


Hydrology 

Surface  hydrology  is  greatly  simplified  in  PRZM  since  it  is  not  the 
main  focus  of  the  model.  Runoff  is  computed  using  the  SCS  curve  number 
method  and  rainfall  that  does  not  run  off,  but  infiltrates.  Evapotranspira- 
tion  is  computed  from  user  inputs  of  pan  evaporation  or  air  temperature. 
Evapotranspiration  is  divided  between  the  plants  and  soil  to  help  track  the 
pesticides.  When  irrigation  is  applied,  infiltration  is  computed  using  the 
Green-Ampt  method.  Percolation  is  based  on  soil  moisture  retention 
capacities.  The  user  may  adjust  the  length  of  time  it  takes  the  soil  to  drain 
to  field  capacity  depending  on  the  soil  texture  and  structure  (Mullins  et  al. 
1996). 
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Pesticide  transport 

PRZM  only  computes  erosion  to  determine  the  amount  of  adsorbed  pes¬ 
ticide  that  leaves  the  field.  Erosion  is  computed  using  the  MUSLE.  Pesti¬ 
cide  movement  considers  volatilization  from  the  soil  as  well  as  from  the 
plant  leaves.  A  balance  between  pesticide  that  is  sorbed  and  in  solution  is 
maintained.  Movement  through  the  vadose  zone  is  controlled  by  Darcy’s 
law  and  Richard’s  equation.  The  transport  equations  can  be  solved  by 
either  finite  difference  or  method  of  characteristics  techniques.  Decay  and 
vapor  phases  of  the  pesticides  are  also  considered  (Mullins  et  al.  1996). 
PRZM  includes  a  Monte  Carlo  mode  that  enables  the  user  to  input  distri¬ 
butions  for  parameters  and  determine  the  probabilities  that  pesticides  will 
leach  below  the  root  zone  (Mullins  et  al.  1996). 


Model  use 


PRZM  has  been  tested  and  validated  in  several  areas  of  the  United 
States  and  is  used  by  EPA  and  chemical  companies  to  evaluate  pesticide 
leaching  (Donigian  and  Huber  1990).  Support,  documentation,  and  PC  ver¬ 
sions  of  the  model  are  available  though  the  EPA’s  CEAM  (Mullins  et  al 
1996). 
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3  Conclusions  and  Model 
Recommendations 


Two  models  are  the  most  comprehensive  of  those  reviewed,  HSPF  and 
SWAT.  The  instream  processes  are  better  in  HSPF,  but  SWAT  deals  more 
completely  with  field  processes  and  groundwater.  Several  areas  could  be 
improved  in  both  SWAT  and  HSPF.  Runoff  in  SWAT  is  currently  based  on 
the  SCS  curve  number  method  and  the  rational  method  to  compute  peak 
runoff  rate.  The  use  of  the  SCS  curve  number  for  computing  water  bal¬ 
ance  is  unacceptable.  This  should  be  improved  by  simulating  breakpoint 
rainfall  and  computing  infiltration  as  in  WEPP.  The  runoff  could  then  be 
routed  using  kinematic  wave  approximations  as  in  HSPF,  WEPP,  or  the 
second  option  in  CREAMS/GLEAMS.  This  could  be  included  as  a  user- 
selected  option  since  use  of  curve  numbers  may  suffice  for  larger  water¬ 
sheds.  HSPF  could  use  updates  to  the  hydrology  section  and  development 
of  an  interface  to  improve  both  preprocessing  and  postprocessing 
capabilities. 

Currently,  sediment  yield  in  SWAT  is  computed  using  the  MUSLE. 

Two  options  exist  for  improving  the  sediment  yield  module  of  SWAT.  The 
first  option  is  to  include  the  option  of  using  the  USLE  or  the  Foster- 
Onstad  USLE  in  addition  to  the  MUSLE.  This  will  allow  the  user  to  select 
the  main  erosion  energy  input.  The  second  option  is  to  include  the  more 
physical  detachment,  transport,  and  deposition  equations  used  in 
CREAMS/GLEAMS.  This  would  increase  computational  time  and  require 
the  hydrologic  modifications  mentioned  previously. 

The  field  and  lake  processes  for  chemical  transport  in  SWAT  are  ade¬ 
quate,  but  little  consideration  is  given  to  channel  processes  other  than 
deposition.  Since  the  time  in  channels  during  transport  is  small  compared 
with  residence  time  in  the  field  or  in  lakes,  this  is  only  a  minor  concern 
for  small  watersheds.  However,  for  large  watersheds  or  basins,  it  repre¬ 
sents  a  problem.  Applying  SWAT  to  an  urban  watershed  may  present  an 
additional  problem.  SWRRB  successfully  simulated  the  effects  of  urbani¬ 
zation  on  a  watershed,  but  at  this  time  SWAT  has  not  been  validated  on  an 
urban  watershed.  Perhaps  an  approach  similar  to  HSPF,  using  different 
modules  to  compute  loadings  from  urban  areas,  could  be  used. 
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^  HSPF  is  selected  as  the  starting  point,  several  upgrades  in  the  hydro- 
logic  portion  of  the  model  are  needed.  The  hydrology  in  HSPF  is  based  on 
the  Stanford  Watershed  Model  developed  in  the  early  1960’s.  Many  of  the 
components  in  use  at  that  time,  although  sufficient,  are  outdated.  Better 
methods  to  compute  evapotranspiration  need  to  be  included;  Penman- 
Monteith,  Hargreaves-Samani,  or  Priestly-Taylor  should  be  considered. 
The  infiltration  section  should  be  upgraded  with  Green-Ampt  methodol¬ 
ogy.  Percolation  and  groundwater  flow  need  to  be  addressed  better.  Appli¬ 
cation  of  nonprecipitation  (irrigation)  events  needs  to  be  included  in  the 
water  balance. 

Sediment  production  simulations  are  adequate,  and  temporal  changes 
in  cover  are  addressed  in  the  AGCHEM  yield  uptake  algorithm.  The 
AGCHEM  module  utilizes  monthly  factors  and  plant/harvest  dates.  Nitro¬ 
gen  fixation  has  been  included  as  well  as  volatilization  and  soluble 
organics.  The  effects  on  groundwater  quality  should  be  included.  How¬ 
ever,  the  inclusion  of  a  groundwater  component  could  be  a  major  undertak¬ 
ing  unless  a  very  simple  approach  is  used.  In  addition  to  these  items,  an 
interface  with  a  GIS  is  needed  as  well  as  a  weather  generator  to  reduce 
data  collection  and  entry  time. 

Current  research  being  undertaken  by  the  EPA  focuses  on  the  develop¬ 
ment  of  a  GIS  and  preprocessor  and  postprocessor  for  HSPF.  The  efforts 
of  the  EPA  have  created  a  software  package  known  as  Better  Assessment 
Science  Integrating  a  Point  and  Nonpoint  Sources  (BASINS).  Current  limi¬ 
tations  of  the  BASINS  package  are  that  many  of  the  HSPF  algorithms  are 
not  available.  Furthermore,  the  database  is  structured  such  that  is  not 
easily  modified  for  site-specific  applications.  An  advantage  of  the 
BASINS  program  is  that  several  options  are  available  for  routing  water 
through  stream  channels.  Other  than  the  HSPF  routing  methods,  QUAL2E 

TOXIROUTE  are  available  within  the  system.  Another  system  devel¬ 
oped  by  the  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  is 
the  Watershed  Modeling  System  (WMS).  This  state-of-the-art  system 
allows  for  utilization  of  available  GIS  data  for  delineation  of  watershed 
boundaries  and  assignment  of  attributes  such  as  land  use  type,  soil  type, 
etc.  The  WMS  offers  excellent  preprocessing  and  postprocessing  capabili¬ 
ties  and  the  flexibility  to  build  unique  watershed  databases  from  user  pro¬ 
vided  data.  The  WMS  also  contains  several  models  including  HEC1, 

TR20,  CASC2D,  and  XMODEL.  BASINS  and  WMS  are  the  two  major 
watershed  modeling  interfaces  that  can  be  used  to  address  total  maximum 
daily  loads,  land  use  changes,  nonpoint  sources  loadings,  etc. 

Upon  reviewing  existing  watershed  water  quality  models,  it  is  recom¬ 
mended  that  the  HSPF  model  be  utilized  as  the  building  block  for  the 
development  of  a  Corps  of  Engineers  modeling  capability  for  predicting 
watershed  runoff  water  quality  as  influenced  by  land  use  change,  develop¬ 
ment,  watershed  improvement  strategies,  and  best  management  practices. 

The  model  should  be  integrated  within  the  existing  WMS  to  facilitate 
preprocessing  and  postprocessing  while  allowing  the  user  maximum  flexi¬ 
bility  to  easily  modify  databases  to  meet  the  demands  of  site-specific 
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applications.  Like  the  BASINS  program,  WMS  is  expanding,  and  it  is 
anticipated  that  other  models  will  be  linked  into  the  system.  CE-QUAL- 
W2  is  listed  among  these  models,  and  this  model  will  offer  an  improved 
reservoir  routing  simulation  option  which  could  be  used  in  conjunction 
with  HSPF.  Another  advantage  of  HSPF  is  that  the  code  is  constantly 
being  updated  and  maintained.  Both  the  EPA  and  the  USGS  provide 
significant  financial  resources  to  ensure  the  code  is  state  of  the  art.  The 
official  version  of  the  code  is  provided  by  the  EPA.  Through  incorpora¬ 
tion  of  the  HSPF  code  into  the  WMS,  the  Corps  of  Engineers  will  have 
the  modeling  capability  required  to  determine  basin-wide  impacts  on  the 
water  quality  of  receiving  waters. 
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